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We report on results of our theoretical study of the in-plane infrared conductivity of the high- 
Tc cuprate superconductors using the model where charged planar quasiparticles are coupled to 
spin fluctuations. The computations include both the renormalization of the quasiparticles and the 
corresponding modification of the current-current vertex function (vertex correction), which ensures 
gauge invariance of the theory and local charge conservation in the system. The incorporation of the 
vertex corrections leads to an increase of the total intraband optical spectral weight (SW) at finite 
frequencies, a SW transfer from far infrared to mid infrared, a significant reduction of the SW of the 
superconducting condensate, and an amplification of characteristic features in the superconducting 
state spectra of the inverse scattering rate 1/r. We also discuss the role of selfconsistency and 
propose a new interpretation of a kink occurring in the experimental low temperature spectra of 
1/r around 1000 cm~^. 
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I. INTRODUCTION 

Most of the existing calculations of the frequency de- 
pendent conductivity a in high-Tc cuprate superconduc- 
tors (for representative examples see Refs. [l|I2H|3HIH3) 
employ the approximation, where interactions of the ex- 
cited states are not taken into account, i.e., the so called 
vertex corrections (VC) in the perturbation expansion for 
a are neglected. It is well known that this approach is not 
gauge invariant and the results may depend on the gauge 
of the vector and scalar potentials. Physically, gauge in- 
variant response is a manifestation of local charge conser- 
vation in the system. We recall that the usual restricted 
conductivity sum rule for the normal state (NS) 
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the so called effective kinetic energy per unit cell (a 
stands for the spin index, is the dispersion relation, 
and rika the occupation factor), is a consequence of the 
gauge invariance-^. For concreteness we have considered 
the conductivity along the x-axis. The NS conductivity 
calculated using a no-VC approach need thus not satisfy 
the sum rule. Similarly, for the superconducting state 
(SCS), the spectral weight of the condensate Ic, which is 
given by 



If+Ic = -{■Ke^/2n^)K , 



(4) 



as obtained within a no-VC calculation, may be rather 
different from that of the corresponding gauge-invariant 
approach. These possible problems call for an assessment 
of the contribution of the relevant VC to If. 

Another reason for studying the VC has been high- 
lighted by Millis and coworkers. If the VC were negligi- 
ble, the infrared (IR) conductivity would be determined 
solely by the quasiparticle self energy, that can be, at least 
in principle, extracted from photoemission data. Millis 
and Drew^ found out, however, that the effective scatter- 
ing rate in optimally doped Bi2Sr2CaCu208+5 obtained 
using the selfenergy estimated from the photoemission 
data, is much higher than the value resulting from the 
IR data. Similarly, Millis et aL^- demonstrated that the 
local selfenergy obtained by fitting (without VC) the IR 
data of thin films of Pr2-a;CexCu04+5 leads to much 
lower (by a factor of 3-4) values of the quasiparticle veloc- 
ity than observed by photoemission. These discrepancies 
appear to imply the presence of large VC. 

A brief account of the VC to the NS conductivity of 
cuprate superconductors has been given by Monthoux 
and Pines in one of the pioneering papers on the nearly 
antiferromagnetic Fermi liquid (NAFL) model^^. For this 
particular model, the VC have been found to cause a 
slight increase (about 20%) of resistivity. Current ver- 
tex renormalization within the conserving fluctuation ex- 
change (FLEX) approachi^ has been systematically in- 
vestigated by Kontani and coworkers (see, e.g., Refs. [isl 
and [ij and references therein). They successfully ex- 
plained the anomalous temperature dependence of the 
Hall conductivity in cuprates and estimated the contri- 
bution of the VC to the normal state conductivity. Ben- 
fatto et ali^ and Aristov and Zeyher— discuss the role 
of the VC in the context of optical response of systems 
possessing the d-density wave ground state, that might 
occur in the pseudogap regime of underdoped cuprates. 
We are not aware, however, of any work incorporating 
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the VC in computations of the optical response in the 
superconducting state. 

The purpose of the present study is to explore the 
role played by the VC, both in the NS and in the SCS 
state, within the spin-fermion model, where charged pla- 
nar quasiparticles are coupled to spin fluctuations (SF). 
We discuss the violation, in the absence of VC, of the NS 
sum rule, the impact of the VC on the dc conductivity, on 
the spectral weight of the superconducting condensate, 
and on the shape of (Ti(w). An important finding is that 
the incorporation of the VC leads to a spectral-weight 
shift from far-infrared (FIR) to mid- infrared (MIR). Fi- 
nally, we address the related changes of sharp structures 
in the inverse scattering rate [1/T](a;). 

The rest of the paper is organized as follows. In Sec. II 
we present the basic equations of the theory, the val- 
ues of the input parameters, and some computational 
details. Since part of the formalism has been already 
detailed in the related previous work by Casek et al^ 
we keep the account short. Section [ni] contains our re- 
sults and discussion. The quasiparticle Green's functions 
have been obtained here using fully selfconsistent Eliash- 
berg theory, whereas an approximate non-selfconsistent 
approach, where the BCS propagator is corrected by the 
coupling to the SF, has been adopted in Ref. 0. In or- 
der to clarify the role of selfconsistency, we first neglect 
the VC and compare the results of the two approaches 
fSec. IIIl!S| . In Sec. IIII B] we then focus on the contribu- 
tion of the vertex corrections. The summary and conclu- 
sions are given in Sec. IIVI 



II. THEORY, INPUT PARAMETERS, AND 
COMPUTATIONAL DETAILS 



to the SF. If we write the right hand side of Eq. ^ 
as a convolution, E = g^xsF * G, then the latter ap- 
proach, called the hybrid approach in the following, yields 
E/i — AfcTi -I- 5^XSF * Go, where Qq stands for the BCS 
propagator with the gap A^. This description is partic- 
ularly suitable for situations, where the coupling to the 
SF is not the only cause of superconductivity and pro- 
vides an additional renormalization of the quasiparticles. 
The magnitude of the superconducting gap can be tuned 
by changing the input amplitude Ag of A^. An impor- 
tant advantage of the hybrid approach is that it avoids 
the necessity of numerical continuation from Matsubara 
frequencies to the real axis. 

Within the framework of the linear response theory, 
the optical conductivity is given by the Kubo formula 



i{uj + iS) 



where Kij is the diamagnetic tensor, 
^ 1 >^ d^Ek 
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In case of a simple cubic or square lattice and the nearest- 
neighbor tight-binding dispersion relation, TrK is pro- 
portional to the band energy (kinetic energy), ^ Efc^fccr, 
per unit cell. For this reason, K is also called the ef- 
fective kinetic energy. Further, Ilij{q,uj) in Eq. ^ is 
the retarded correlation function of the paramagnetic- 
current-density operator 

i r°° 

n,,{q,u;) = -J die'"*([j-(q,i),Ji(-q,O)])0(i) . (8) 



Within the framework of the spin-fermion model, the 
electronic quasiparticles are renormalized via a coupling 
to SF. These are described by the spin susceptibility 
XSF{k,E) of the form motivated by neutron data and 
the electronic selfenergy (2x2 matrix) of the most so- 
phisticated version of the theory is given hy^ 
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I](fc,iS) = ^ XSF{k-k'AE -iE')g{k',iE') 



(5) 

Here g is the coupling constant, xsf the Matsubara coun- 
terpart of the spin susceptibility, and Q is the Nambu 
propagator of the renormalized electronic quasiparticles, 
g = [iEtq — {ek — fj,)T3 — 'E,{k,iE)]~^ . The sum runs 
through the Bloch vectors k' in the first Brillouin zone 
{N in total) and the fermionic Matsubara energies i£" = 
i7r(2n -|- !)//?. Note that equation ([5]) is a compressed 
form of the Eliashberg equations. 

A popular approximation to the fully selfconsistent 
treatment of Eq. (O, used, e.g., in Refs. ITHgIITsI consists 
in starting with the BCS Nambu propagator with an es- 
timated superconducting gap of d^i_y2 symmetry and 
including the lowest order correction due to the coupling 



In the lowest order approximation, the correlator cor- 
responds to a simple bubble diagram containing two in- 
dependent quasiparticle lines. In the limit of small q, its 
Matsubara counterpart can be written as 
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n,, (q,iH = -72T7«E _ _ 

xTT[g{k + q,iE + ihiy)g{k,iE)] . (9) 



This is the approximation mentioned in the introduction, 
that is not gauge invariant and may lead to a violation 
of the sum rule H]). 

A general field theoretical method to overcome the 
problem of gauge invariance was constructed by Nambi*i^ 
in the context of the BCS theory. In the present context 
it leads to the replacement of formula ^ with 
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The bare vertex Tj = {de / dkj)To has been replaced with 
a renormahzed one, rj{k,iE,q,ihi') (2x2 matrix). For 
the theory to be gauge invariant, the renormalization of 
the vertex representing the coupHng of the quasiparticles 
to the electromagnetic field, has to be consistent with the 
renormalization of the quasiparticles by the interaction. 
It is known that F has to obey the generalized Ward 
identity, which is in fact a reformulated charge conserva- 
tion law^— . The minimal vertex satisfying the generalized 
Ward identity in the present case of the selfenergy given 
by Eq. ([5]) is the solution to the Bethe-Salpeter equation 

^ ^ k',iE' 

X g{k' + q,iE' + ihiy)Tj{k',iE',q,ihiy)g{k',iE') . 

(11) 

The second term on the right hand side of the equa- 
tion is the VC. The correlator resulting from this gauge- 
invariant approach corresponds to the sum of all ladder 
diagrams, where non-crossing SF lines connecting the two 
quasiparticle lines are inserted into the conductivity bub- 
ble. 

Formally, the hybrid approach can be made gauge in- 
variant as well by using a somewhat simpler form of the 
renormahzed vertex. The results, however, are qualita- 
tively different from those of the selfconsistent approach. 

The input quantities of the theory are the dispersion 
relation e^, the chemical potential (or the electron den- 
sity), the spin susceptibility xsf, and the coupling con- 
stant g. We have used the second nearest neighbor tight- 
binding dispersion relation and the model spin suscepti- 
bility of the same form as in Ref. d, containing the reso- 
nance mode and a continuum with dimensionless spectral 
weights of 0.016a/ and 0.016c, respectively. The values 
of all input parameters are given in Table ID They are 
the same as in Ref. H, except for g (1.52 eV in Ref. @). 
The present value of g of 2 eV yields Tc = 64 K and 
A = 20 meV (A is the amplitude of the gap). In the 
computations of Sec. Ill B, slightly different values of t, 
n, and g have been used, in Table [J they are given in the 
brackets.^l. For g = 2eV (5 = 3eV) we obtain = 77K 
and A = 22 meV (T^ = 89 K and A = 27 meV). 

Note, that the present formulation assumes a single 
band only. In order to compare our results with experi- 
mental data on bilayer cuprates, we multiply the conduc- 
tivity and related quantities (like K) by a factor Np\ = 2 
reflecting the two planes in the unit cell. To address this 
issue carefully, we have also performed calculations with 
a model containing two bands - bonding and antibond- 
ing. The spin susceptibility was redistributed between 
the odd- and even-channels, the resonant mode being ac- 
tive in the odd channel?^. The single-particle part of our 
computations is similar to that of Ref. [2^. We have ne- 
glected the VC for simplicity. The results indicate, that 
for reasonable values of the intrabilayer hopping matrix 
elements (of the order of lOOmeV), the band splitting has 
a negligible effect on the in-plane conductivity, changing 



only the relative contributions of the bonding and the 
antibonding band to the resulting spectra. 

The iterative solution of the selfenergy equation ([5]) 
and especially of the Bethe-Salpeter equation pT|) in- 
volves many convolutions both in Matsubara frequencies 
and in fc-space. Together with the requirement of high 
accuracy data for the subsequent analytical continuation 
to the real axis, this leads to a computationally very de- 
manding task. We have performed the convolutions using 
the FFT algorithm taking the full advantage of the sym- 
metries of E and F. Typically, we have used a grid of 
96 X 96 points in the Brillouin zone and a cutoff of 10 eV 
(approximately 3.5 times the bandwidth) to limit the 
number of Matsubara frequencies. We have checked, by 
varying the density of the grid and the cutoff, that these 
values are sufficient. The analytical continuations to the 
real axis were performed using the standard method of 
Fade approximants— . 



III. RESULTS AND DISCUSSION 
A. Role of selfconsistency 

Here we compare the conductivity spectra computed 
using the (nonselfconsistent) hybrid approach with those 
of the selfconsistent Eliashberg theory. Figure [T] shows 
the spectra of cti for the NS at T = 100 K (a) and for 
the SCS at T = 20 K (b). For the NS, the conductivity 
profiles are quite featureless. The dc value of the hybrid 
approach is by about 30% smaller. This is because the 
hybrid approach, which is - only for the NS - equiva- 
lent to the first iteration of Eq. ([5]), overestimates the 
magnitude of the quasiparticle selfenergy. 

Next we discuss the spectra for the SCS. They display 
an increase of cti starting at low frequencies, which be- 
comes steeper around Tiujq -I- A as discussed in Ref. H. 
The character of the maximum following the onset, how- 
ever, depends on the level of theory used. For the hybrid 
approach, the conductivity exhibits a very broad max- 
imum, whereas the selfconsistent approach yields a rel- 
atively sharp maximum near 800 cm~^, consistent with 
experimental dat a^^i^^'^^ . The origin of this difference 
can be traced back to that of the quasiparticle spectral 
functions, using the method presented in Ref. 6. It is 
based on the formula valid for the NS {q — 0) 
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A{k,E)A{k,E + huj)[nF{E) - npiE - 



huj)]dE , 

(12) 



where A{k, E) is the spectral function and np the Fermi 
function, and a similar formula involving the matrix spec- 
tral function valid for the SCS. These formulas allow one 
to understand the structures in the spectra of ai in terms 
of "transitions" between various components of A. It has 
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a [A] d [A] 


t[cV] 


t'[eV] n 


5[eV] 


Tiloq [eV] 


r[evi CM 


RLJc[eV] rc[eV] CcM 6m 6c 


3.828 11.650 


0.250 (0.350) 


-0.100 0.76(0.82) 


2.0(2.0/3.0) 


0.040 


0.010 2.35 


0.400 1.000 0.52 1 4 



TABLE L Values of the parameters used in the computations. Lattice parameters are denoted by a and d, t and t' are the 
parameters of the dispersion relation, n is the number of electrons per unit cell, the meaning of the parameters uoo, T, C, ijJc, 
Vc, and Cc specifying the spin susceptibility is explained in Ref. 0. For the normal state, we take F = 0.070 eV and C ~ 1.57a. 
The value of the input amplitude Ao of Afc used within the hybrid-approach was adjusted to yield the same amplitude of the 
resulting gap A as the selfconsistent approach. In the computations of Sec. Ill B, the values in the brackets have been used, 
the spectra shown in the figures E] HI O |6] have been obtained with g — 3 eV. 
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FIG. 1: Real part of the optical conductivity as obtained by 
using the selfconsistent approach (solid lines) and the hybrid 
approach (dashed lines) for the normal state at T = 100 K 
and for the superconducting state at T = 20 K. 
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FIG. 2: Quasiparticle spectral functions for selected fc-points 
defined in Ref. 12 and shown in the insets as obtained by 
using the selfconsistent approach (solid lines) and the hy- 
brid approach (dashed lines) for the superconducting state at 
T = 20 K. Our notation of spectral structures is introduced 
in panels (c) and (d). 




been shown in Refs. that the onset of ai in the SCS 
including the maximum is determined by transitions be- 
tween the quasiparticle peaks and the incoherent parts 
of A. The spectral functions for selected fc-points are 
presented in Fig. [2l A detailed discussion of the spectral 
structures obtained at the hybrid level can be found in 
Ref. 0. Here we concentrate on the comparison between 
the results of the selfconsistent and the hybrid approach. 
The positions of the quasiparticle peaks Q, Q' are identi- 
cal, which has been achieved by tuning the input BSC gap 
of the hybrid approach. The maxima Pi and P[ of the in- 
coherent part in the selfconsistent case are located closer 
to the quasiparticle peaks and are considerably sharper. 
This causes the steeper onset and the sharper maximum 
of (Ti. It can be seen that this maximum simply reflects 
the maxima of the incoherent part of A - Pi and P[. 



B. Vertex corrections 

In this section we focus on the central topic of the 
paper of how the model spectra change when the VC 
are included. Figure [3] shows the optical conductivity 
and three related quantities: the inverse scattering rate 
[1/t\{lo) and the mass enhancement factor [m* /m\{uj)^ 
that are defined by the extended-Drude-model formula 

^'^^^ = h / u ^ — ■ r */ u ^ ' (^^^ 
[1/t\(lo) — \ijj[m* I m\(Lo) 

where Wpi is the plasma frequency, and W(uj) — 
(27r)-M2|a;/r(w)]/dw2. For the NS of a weakly cou- 
pled isotropic electron-phonon system the function W{uj) 
is approximately equal to the electron-phonon spectral 
density^*'^^. The solid and the dashed lines correspond 
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FIG. 3: The spectra of the optical conductivity, the inverse 
scattering rate, the mass enhancement factor, and the func- 
tion W calculated with the vertex corrections neglected (NV) 
and with the vertex corrections included (VC) for the set of 
input parameters leading to realistic values of Tc and A. The 
insets of panel (a) shows the vertically zoomed area near the 
conductivity maximum and the spectra of (T2. In panel (d), 
only the superconducting state spectra are shown. The ver- 
tical dashed lines are drawn to guide the eye. A notation of 
spectral structures is introduced in (b). 



to the SCS at T = 20K, the dashed-dotted and the dotted 
lines to the NS at T = 100 K. The sohd and the dashed- 
dotted lines (the dashed and the dotted lines) correspond 
to the selfconsistent computation with (without) the VC. 
The values of the dc conductivity CTdc, the effective kinetic 
energy K (multipUed by d to obtain the dimension of en- 
ergy), the spectral weight at finite frequencies //, and 
that of the condensate, /c, are given in Table |TT1 The 
values of // and 1^ have been multiplied by 2h?d/'Ke^ to 
allow a direct comparison with those of Kd. 

We begin our discussion with the NS. It can be seen in 
Table HIl that in the absence of VC the sum rule If ^ K 
(III) is not fulfilled. It means that the conductivity pos- 
sesses an unphysical singular component with the spec- 
tral weight Ic determined by Eq. (|4]). The values of Ic 
of ca 3% (ca 5%) of -K for g = 2 eV (g = 3 eV) are 
small but significant. The unphysical component mani- 
fests itself also in the spectra of related quantities - see 
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If [eV] Ic/\K\ [%] adc[n" 
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NV 


0.27079 0.17609 


0.43205 3.2 1.62 


• 10^ 


VC 


0.29752 0.14936 


0.44624 0.1 1.41 


■ 10* 


g = 


3 eV = 89 K 


Ax = 27 meV (20 K) 






T = 20K 


T = 100 K 






Kd = -0.40741 eV 


Kd = -0.40699 eV 






If [eV] h [eV] 


If [eV] IJ\K\ [%] ai,[Q- 


^cm-i] 


NV 


0.28871 0.11870 


0.38580 5.2 8.0 


10^ 


VC 


0.31312 0.09429 


0.40699 0.0 7.5 


10^ 



TABLE II: Values of the effective kinetic energy K, the spec- 
tral weight at finite frequencies //, and the spectral weight 
of the singular component (condensate) Ic computed with- 
out vertex corrections (NV) and with the vertex corrections 
included (VC). The data are presented for two values of the 
coupling constant g, for the superconducting state at T = 20K 
and for the normal state at T = 100 K. For the normal state, 
the values of the dc conductivity are also given. 



the spectra of (72 in the inset of Fig. [3] (a), the drop of 
1/t at low frequencies [Fig. [21(b)], and the corresponding 
divergence of m* /m [Fig. [3] (c)]. With the VC included 
the sum rule is satisfied. The spectral weight increase 
due to the VC [A//(VC)] is equal to Ic (within the nu- 
merical error related to discrete fc-sampling). Based on 
the relatively small values of Ic, the changes of the spec- 
tra due to the VC can be expected to be small. Indeed, 
only a slight decrease of cti in the FIR can be observed 
in Fig. [21 The dc conductivity also decreases (see Table 
HI)) which is consistent with the results of Monthoux and 
Pincs^^. This trend, however, is not universal, as will be 
discussed in Sec. IIII B II 

In the SCS, the VC increase //, which leads to a re- 
duction of Ic (see Table [ill • This effect is explored in 
detail in Sec. IIII B 21 The real part of the conductivity 
is affected mainly near then maximum around 800 cm^^: 
the VC make it slightly sharper. The changes are further 
amplified in the spectra of related quantities 1 /r and W 
to be discussed in Sec. IIII B 31 



1. Spectral weight transfer from FIR to MIR 

Figure ^a) shows a decrease of fii in FIR with the 
incorporation of the VC which seems to be inconsistent 
with the increase of the finite-frequency spectral weight 
//. A resolution of this apparent controversy is provided 
by Fig. [H which shows the contribution Acti (VC) of the 
VC to CTi in a wide spectral range. It can be seen that 
A(Ji(VC) is negative for ut < 2000 cm^^ but positive for 
u! > 2000 cm~^. The magnitude of the contribution of 
the latter spectral range to // is larger than that of the 
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FIG. 4: (a) Contribution of the VC to the real part of the 
normal state conductivity at T = 100 K for various values of 
the relative spin-continuum weight be- The low energy parts 
of the spectra are shown in the inset, (b) The same for the 
superconducting state at T = 20 K. 
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FIG. 5: (a) The coupling-constant dependence of the effective 
kinetic energy and of the intraband spectral weight obtained 
with the vertex corrections neglected (NV) and with the ver- 
tex corrections included (VC). (b) The coupling-constant de- 
pendence of the ratio of the condensate weight obtained with 
the vertex corrections included to that obtained using the 
vertex corrections neglected. 



former, which leads to the total increase of //, both for 
the NS and for the SCS. For the SCS, the magnitude 
of Ao'i(VC) becomes small as the frequency approaches 
zero. This is due to the vanishing real part of conductiv- 
ity in this limit. Figure [4] further demonstrates that the 
maximum of A(Ti(VC) shifts towards higher frequencies 
and becomes broader with increasing weight of the spin 
fluctuation continuum. 

In order to explain the behavior displayed in Fig. |4l 
we use an extension of the formula (jl2p that ap- 
plies to the theory involving the VC. In the general- 
ized formula, which has been obtained by manipula- 
tions starting from Eq. pO]) . the integral on the right 
hand side of Eq. (fT^ is replaced with a more gen- 
eral convolution of the spectral functions of the form 
/ dE' J dE"A{k, E')A{k, E")E{E', E", huj), with the the 
Matsubara counterpart of the convolution kernel given by 

EiE',E",in.) = -ImV , miEAh.) ^ ^ 

^ ' ' hp ^ {iE + \hv - E'MiE - E") 

(14) 

If we use the bare vertex F — de/dk, we arrive at the re- 
lation S (X S{E"+huj-E')[nFiE")-nF{E" + hLj)] which, 
after performing one of the integrations in the convolu- 
tion, provides exactly the formula (fT^ . In contrast, any 
kernel S corresponding to a frequency dependent renor- 
malized vertex will be broader than the delta function 
and will result in a broadening of the conductivity pro- 
file. 

The effect of the VC on the conductivity can be vaguely 
viewed as consisting of two ingredients: (a) the increase 
of the overall spectral weight // and (b) the broaden- 
ing of the conductivity profile, which causes a transfer 
of spectral weight towards higher frequencies. The point 
(b) accounts for the trends shown in Fig. 21 The change 



of the dc conductivity is determined by a competition of 
(a) and (b). For our choice of the input parameters, the 
VC reduce the dc conductivity, in some cases, however, 
where A7y(VC) is large, an increase of iTdc may occur. 
Some examples can be found in Ref. 

The spectral weight redistribution of Fig. 0] is similar 
to but considerably smaller than the one that could be 
expected based on the arguments by Millis et a^. Note, 
however, that the discrepancy between the value of the 
quasiparticle velocity along the Brillouin zone diagonal 
resulting from the present computations of 1.1 eVA for 
g = 2 eV (0.75 eVA for 5 = 3 eV) and the experimental 
value for YBCO of ca 1.6 eVA^ is not as dramatic as in 
Ref. M- 



2. Effect of the vertex corrections on the spectral weight of 
the condensate 



We have already noted that the VC significantly in- 
crease the values of // and reduce those of Ic- This is 
further documented in Fig. [5l which shows // and Ic as 
functions of the coupling constant g. For values of g^ of 
4—10 eV^, leading to reasonable values of Tc and A, the 
magnitude of A//(VC) is about 6% of that of K. It can 
be seen in part (b) that the corresponding change in the 
condensate weight Ic is up to 20%. At low values of g, 
A//(VC) scales with g^. This can be easily interpreted 
by considering the diagram series for the correlator (|10p . 
The contribution of a diagram with A'^ spin-fluctuation 
lines connecting the quasiparticle lines in the bubble is 
proportional to g^^. In the limit of small g, only the 
lowest-order diagrams {N — and = 1) survive, lead- 
ing to the observed behavior. 
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The temperature dependence of the effective kinetic 
energy K and the band energy KE = ^ ^k^ka ob- 
tained within EUashberg theory with SF, was aheady 
discussed in Ref. [sil . Our approach yields a similar be- 
havior of the two quantities. The resulting value of 
K{m,T = 20 K) - K{SCS,T = 20 K) is also similar 
to that obtained by Casek et al.^ using the hybrid ap- 
proach. The analysis published in Ref. 6^, however, can- 
not be easily extended to the present fuUy-selfconsistent 
theory. The contribution A//(VC) of the VC to // is 
only weakly temperature dependent above T^. Below 
Tc, A//(VC) slightly decreases, exhibits a minimum, and 
then increases. The low temperature value is somewhat 
higher than that of the NS. 



3. Spectral structures of [1/t\(ijj) and W{uj) for the 
superconducting state 

Having discussed the global trends of the spectral 
changes due to the VC, we concentrate here on pro- 
nounced features in the spectra of l/r and W. The SCS 
spectra of l/r shown in Fig. [3] exhibit the famous onset 
starting around the frequency of the resonance mode wqi 
becoming steeper around 500cm^^ (this feature is labeled 
as 5), and reaching a sharp maximum around 800 cm^^. 
Surprisingly, the maximum is followed by a kink labeled 
as K. Note that the spectra are fairly similar to the ex- 
perimental ones of optimally doped material o^^i^^'^^ . All 
these features appear already in the spectra of cti, how- 
ever, they are more pronounced in those of l/r. The 
function W is approximately proportional to the second 
derivative of l/r and can thus be expected to possess two 
maxima corresponding to the structures S and K and a 
minimum close to the maximum of l/r. This is indeed 
the case, as shown in Fig. [3] (d). 

The origin of the feature S has been elucidated by 
Casek et . It is due to the appearance above Tiloq -f A 
of excitations of the nodal region, consisting of a nodal 
quasiparticle, an antinodal quasiparticle, and the reso- 
nance mode. The arguments of Ref. even though for- 
mulated at the level of the hybrid approach, remain to 
be valid also in the context of the present fully selfcon- 
sistent theory. As noted for the first time by Carbotte, 
Schachinger, and Basov^, the energy of the relevant max- 
imum in W is close to Tlujq + A. 

The sharp maximum of l/r and the corresponding 
minimum of W can be shown to result from transitions 
Pi Q' and Q ^ P[ (see Fig. H]). The energy of the 
structure is close to Tiujq + 2 A, which can be understood 
using the arguments presented in Ref. [^. The presence of 
this characteristic energy scale in the optical spectra of 
the high- Tc cuprates has been for the first time predicted 
by Abanov, Chubukov, and Schmaliark^.. In their work, 
the corresponding spectral structure is attributed to pro- 
cesses involving a Bogoljubov quasiparticle with energy 
A and a sharp onset of the incoherent part of A{k, E) at 
A + fujjQ. 
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FIG. 6: The function W for the superconducting state at 
T = 20K obtained with the vertex corrections neglected (NV) 
and with the vertex corrections included (VC) for three values 
of the energy of the resonant mode in the spin susceptibility. 
The vertical lines are located at the characteristic energies of 
hujo + A, h,u)o -f 2A, and 2S,ujo -I- 2A. 

The kink labeled as K can also be interpreted in terms 
of the quasiparticle spectral functions. Note that the 
maxima Pi and P[ in Fig. [5] that are connected to the 
maximum of di and the related structures in the spectra 
of l/r and VF, are followed by shoulder features, labeled 
as P2 and on their high-energy sides. They can be ex- 
pected, based on Eq. to manifest themselves also in 
the conductivity. Our detailed calculations show that (a) 
the shoulder features are indeed responsible for the kink 
K \ii l/r and (b) they can be attributed to excited states 
involving two magnetic excitations. We recall that the 
maxima Pi and P[ correspond to excited states involv- 
ing one Bogoljubov quasiparticle and just one magnetic 
excitation. The characteristic energy of the shoulder fea- 
tures is 2T1UJ0 + A and that of the kink K (which can 
be associated with transitions P2 Q' and Q ^ P2) 
is 2fiu)Q -\- 2A. Abanov, Chubukov, and SchmalianS also 
predicted a structure located at 2TujJq + 2A. Their in- 
terpretation, however, differs from ours. In their theory, 
the structure is due to transitions between negative- and 
positive- energy satellites (incoherent components) of A. 

In order to check the proposed assignment of the spec- 
tral structures, we have studied the wo-dependence of l/r 
and W . Figure [6] shows the spectra of W for three values 
of the energy of the magnetic mode. It can be seen that 
the first maximum and the minimum are located close to 
TiLjQ + A and Tilo^ + 2A, respectively, in agreement with 
the above considerations. The second maximum corre- 
sponding to the kink is, in the absence of VC, located 
somewhat below 2TujJq -|- 2A. 

Finally, we address the role of the VC. They lead (a) 
to an increase in the amplitude of the structures of l/r 
and VF, which is due to the combined effect of the sig- 
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nificant decrease of (Ji [see the inset of Fig. [3] (a)] and 
a sharpening of the features in cti, and (b), to a shght 
shift of the structures towards lower energies. The shift 
may be partiaUy caused by a weak attractive interaction 
between the quasiparticles leading to an excitonic effect. 

IV. SUMMARY AND CONCLUSIONS 

The changes of the infrared conductivity caused by the 
vertex corrections (VC) are not dramatic, which provides 
a justification for earlier computations, where these cor- 
rections were neglected. Some aspects of the changes, 
however, appear to be important. 

(i) The normal state conductivity computed without the 
VC does not satisfy the restricted sum rule. The calcu- 
lated spectral weight at finite frequencies (//) is by a few 
percent lower than the value dictated by the sum rule. 
The incorporation of the VC leads to the increase of // , 
removing this discrepancy. 

(ii) The increase of // is associated with a broadening 
of a\. The far-infrared conductivity (including the dc 
value) may decrease or increase depending on the magni- 
tude of the change of // and the degree of the broadening; 
at high frequencies the conductivity increases. For our 
values of the input parameters, a\ decreases (increases) 
below (above) ca 2000 cm~^. 

(iii) The increase of // occurs also for the superconduct- 
ing state. Since the sum of /j and the spectral weight of 
the superconducting condensate remains constant, this 
implies a reduction of the latter. For relevant values of 
the input parameters, the weight of the condensate de- 
creases by 15-20%. 

(iv) The VC lead to an amplification of the characteristic 
features in the superconducting state spectra of the in- 
verse scattering rate 1/t, that have been used to support 
the spin-fluctuation scenario: the onset of 1/r around 



500 cm~^ becomes steeper and the maximum around 
800 cm~^ more pronounced. 

In addition to studying the changes brought about by 
the VC, we have also investigated the role of self consis- 
tency by comparing the results obtained using the self- 
consistent Eliashberg equations with those of the nonself- 
consistent hybrid approach (i.e., approximately, the first 
iteration of the equations). The main results are (a) the 
hybrid approach considerably overestimates the magni- 
tude of the quasiparticle selfenergy, which leads to lower 
values of the conductivity in far-infrared, and (b) some 
spectral features, in particular the sharp maximum in the 
superconducting state spectra centered at ca 800 cm~^ 
and the kink at ca 1000 cm~^ appear only at the selfcon- 
sistent level. With the aid of the quasiparticle spectral 
function A, we attribute the kink to the onset of tran- 
sitions involving incoherent satellites of A corresponding 
to states with doubly excited resonance mode.'^S, 

The computed spectra are in reasonable agreement 
with experimental data of optimally doped cuprates, in- 
cluding such details as the shape of the maximum in the 
inverse scattering rate. In addition, the theory allows one 
to interpret most of the features of the data in terms of 
Bogoljubov quasiparticles and magnetic excitations. In 
the superconducting state, a crucial role is played by the 
resonance mode. We are not aware of any comparable 
interpretation in terms of electron-phonon coupling. 
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